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Sharp  quantum  phase  transitions  typically  require  a  large  system  with  many  particles.  Here  we 
show  that  for  a  frustrated  fully-connected  Ising  spin  network  represented  by  trapped  atomic  ions, 
the  competition  between  different  spin  orders  leads  to  rich  phase  transitions  whose  sharpness  scales 
exponentially  with  the  number  of  spins.  This  unusual  finite-size  scaling  behavior  opens  up  the 
possibility  of  observing  sharp  quantum  phase  transitions  in  a  system  of  just  a  few  trapped  ion  spins. 

PACS  numbers: 


Quantum  simulators  are  motivated  by  the  promise  of 
gaining  insight  into  many-body  quantum  systems  such 
as  high-Tc  superconductors  or  complex  arrangements  of 
interacting  spins.  Cold  atomic  systems  form  a  promising 
platform  for  quantum  simulation,  as  the  interactions  be¬ 
tween  particles  can  be  under  great  control.  A  good  exam¬ 
ple  is  a  collection  of  trapped  and  laser-cooled  atomic  ions, 
each  representing  an  effective  spin  that  can  be  made  to 
interact  with  all  the  others  by  modulating  the  Coulomb 
interaction  between  ions.  By  applying  spin-dependent 
optical  dipole  forces  it  has  been  shown  that  a  crystal  of 
ions  provides  an  ideal  platform  to  simulate  intractable 
interacting  spin  models  0,1-  Following  this  proposal, 
recent  experiments  have  simulated  quantum  magnetism 
with  a  few  ions  M.  For  three  or  more  ions,  the  long- 
range  coupling  between  the  spins  can  provide  frustrated 
interaction  patterns  or  competition  between  various  spin 
orders  and  interesting  quantum  phase  diagrams  p,  6j. 
However,  the  observation  of  a  quantum  phase  transition 
typically  requires  a  large  system  with  many  particles,  as 
the  width  (sharpness)  of  a  quantum  phase  transition  usu¬ 
ally  scales  with  1  /TV,  the  inverse  of  the  number  of  par¬ 
ticles  [?].  With  such  slow  finite-size  scaling  laws,  small 
networks  of  trapped  ions  realized  in  current  experiments 
(N  <  20)  are  not  expected  to  exhibit  sharp  transitions 
between  distinct  quantum  phases. 

In  this  paper,  we  show  the  surprising  result  that  sharp 
phase  transitions  can  indeed  be  observed  with  just  a  few 
atomic  ions.  This  is  due  to  unusual  finite  size  scaling 
laws  in  this  frustrated  spin  network,  where  the  sharpness 
of  some  phase  transitions  scales  exponentially  instead  of 
linearly  with  1/N.  By  controlling  a  single  experimen¬ 
tal  parameter  that  determines  the  pattern  of  spin-spin 
couplings  between  the  ions,  we  show  that  the  expected 
ground  state  emerges  from  a  delicate  compromise  be¬ 
tween  the  couplings.  Frustration  in  the  spin  network 
leads  to  a  variety  of  spin  orders,  with  the  number  of  dis¬ 
tinct  phases  increasing  rapidly  with  the  number  of  ions. 
We  construct  the  complete  phase  diagram  for  small  spin 
networks  realizable  with  the  current  technology.  The 
sharp  phase  transition  is  characterized  in  detail  with  an 
explanation  of  its  unusual  finite-size  scaling  behavior. 

We  consider  a  small  crystal  of  ions  confined  in  a 


one-dimensional  harmonic  trap.  The  spin  states  of  the 
ions  are  represented  by  two  internal  states,  referred  as 
|t)  and  |t),  and  the  effective  spin-spin  interaction  be¬ 
tween  the  ions  is  induced  with  off-resonant  bichromatic 
laser  beams  Hi!.  The  ion-laser  coupling  Hamilto¬ 
nian,  written  in  the  rotating  frame,  has  the  form  H  = 
[^cos (Skxn  +  /rt)cr^  +  0?  where  Q  is  a  Ra¬ 

man  Rabi  frequency,  5k  is  the  wave  vector  difference  be¬ 
tween  the  two  Raman  beams  (which  is  assumed  to  be 
along  the  radial  direction  x) ,  fi  is  the  beatnote  or  detun¬ 
ing  between  the  two  laser  beams,  and  a %  are  Pauli  ma¬ 
trices  describing  the  spin  of  the  nth  ion,  and  B  is  an  effec¬ 
tive  magnetic  field  induced  by  radiation  that  coherently 
flips  the  spins.  In  the  rotating  frame,  the  radial  coordi¬ 
nate  xn  is  expanded  in  terms  of  the  transverse  phonon 
modes  ak  as  xn  =  J2k  k/  {2mujk){a\eluJkt  +  ake~luJkt), 
where  m  is  the  atomic  mass,  uvk  is  the  eigenfrequency  of 
the  kth  normal  mode  of  the  ion  crystal,  and  is  the 
eigenmode  transformation  matrix.  We  use  transverse 
phonon  modes  because  they  can  more  easily  be  scaled 
up  to  large  systems  1,1.  Under  the  Lamb-Dicke  crite¬ 
rion  r]n:k  =  b^dk^h/ (2mojk)  1,  the  Hamiltonian  H 
is  simplified  to  H  =  -HUt  J2nk  On,k  sin  {pit)  a*(alez“kt  + 

If  we  assume  that  the  laser  detuning  /a  is  not  resonant 
with  any  phonon  mode  with  the  condition  \u)k  —  n\  ^ 
T]n,k^  satisfied  for  all  n  modes  k ,  the  probability  of  ex¬ 
citing  any  phonon  mode  |U^n?/c6^/2  ( uvk  —  /i)|2  is  negligi¬ 
ble.  We  can  therefore  adiabatically  eliminate  the  phonon 
modes  and  arrive  at  the  following  effective  spin-spin  cou¬ 
pling  Hamiltonian  0® 

Hs=Y,JmnVZm<  +  Bj2<’  W 

m,n  n 

where  the  coefficients 

_  (hQ5k)2 
mn“  2m 

This  Ising  Hamiltonian  is  a  pillar  of  many-body 
physics,  and  its  properties  have  been  exhaustively  stud¬ 
ied  under  various  conditions  [7].  For  instance,  the  ground 
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Figure  1:  (Color  online)  Illustration  of  spin  orders  in  two  frus¬ 
trated  Ising  networks  with  competing  long-range  interaction 
for  N  =  7  ions.  Left  panel :  Ferromagnetic  (FM)  order  with 
detuning  jl  —  5.1;  Right  panel :  Kink  order  with  jl  =  5.3.  The 
thickness  of  the  edge  in  the  graph  represents  the  strength  of 
each  coupling.  Positive  (AFM)  spin  couplings  are  indicated 
in  red,  while  negative  (FM)  couplings  are  indicated  in  black. 


state  of  the  Ising  Hamiltonian  is  well  understood  when 
the  coupling  coefficients  Jmn  are  uniform,  or  nonzero 
only  for  nearest  neighbors.  However,  here  we  have  an  ex¬ 
tended  Ising  network  where  the  coupling  coefficients  Jmn 
are  inhomogeneous  (both  in  magnitude  and  sign)  and  ex¬ 
tend  over  long  range  m-  The  strong  competition  among 
these  interaction  terms  (even  with  B  =  0)  will  generally 
lead  to  highly  frustrated  ground  states  where  individual 
bonds  are  compromised  in  order  to  reach  a  global  en¬ 
ergy  minimum.  For  arbitrary  coupling  coefficients  Jmn, 
the  Hamiltonian  (1)  generally  belongs  to  the  complexity 
class  of  NP-complete  problems  0,  meaning  that  calcu¬ 
lating  attributes  such  as  the  system  ground  state  becomes 
intractable  when  the  system  size  is  scaled  up. 

We  consider  the  case  where  the  coupling  coefficients 
Jmn  are  controlled  by  a  single  experimental  parameter, 
the  laser  detuning  /i  [J-Q.  To  determine  Jmn  from  any 
detuning  q  with  the  formula  (2),  we  first  need  the  normal 
mode  transformation  matrix  b^.  This  is  obtained  by  find¬ 
ing  the  equilibrium  positions  for  a  given  number  of  ions  in 
a  harmonic  trap  and  then  diagonalizing  the  Coulomb  in¬ 
teraction  Hamiltonian  expanded  about  the  ions’  equilib¬ 
rium  positions.  With  a  single  control  parameter  /q  we  are 
not  able  to  program  arbitrary  coupling  coefficients  Jmn. 
However,  the  interaction  pattern  is  sufficiently  complex 
to  allow  frustrated  ground  state  configurations  and  rich 
phase  transitions.  To  illustrate  this  we  show  in  Fig.  1  a 
coupling  pattern  for  N  =  7  ions  and  its  associated  ground 
state  spin  configuration  for  B  =  0.  The  coupling  pattern 
is  represented  by  a  graph  where  the  color  and  the  thick¬ 
ness  of  each  edge  represents  respectively  the  sign  (ferro¬ 
magnetic  or  antiferromagnetic)  and  the  magnitude  of  the 
coupling.  In  Fig.  1(a),  we  find  a  ferromagnetically  or¬ 
dered  ground  state  with  all  the  spins  pointing  to  the  same 
direction.  However,  in  this  ferromagnetic  state,  some  of 
the  bonds,  such  as  the  strong  antiferromagnetic  bond  be¬ 
tween  the  ions  1  and  7,  are  compromised,  and  due  to  this 
frustration,  the  ground-state  spin  configuration  is  very 
sensitive  to  the  strength  of  the  coupling.  If  we  adjust  the 
detuning  g  by  a  small  fraction  of  the  trap  frequency,  the 
ferromagnetic  bonds  of  the  ion  pairs  (1,  5)  and  (3,  7)  are 


slightly  weakened  (see  Fig.  1(b))  and  the  antiferromag¬ 
netic  bond  (1,7)  dominates  and  flips  the  spin  direction 
of  the  entire  left  (or  right)  half  of  the  ion  crystal.  This  is 
a  phase  transition  from  ferromagnetic  order  to  a  “kink" 
order,  with  a  kink  in  the  spin  direction  between  the  4th 
and  5th  ions  counting  from  either  the  left  or  the  right 
side. 

To  show  the  rich  phase  diagram  for  this  system,  in  Fig. 
2(a)  we  list  all  different  spin  phases  at  B  =  0  for  a  small 
Ising  network  with  3,  5,  7  and  9  ions.  For  an  odd  number 
of  ions,  the  phase  diagram  is  more  interesting  and  fea¬ 
tures  a  larger  variety  of  spin  orders,  because  the  left-right 
reflection  symmetry  in  a  linear  ion  crystal  can  be  sponta¬ 
neously  broken.  In  Fig.  2,  for  convenience,  we  denote  the 
laser  detuning  q  with  a  re-scaled  dimensionless  parameter 
q  by  labeling  the  phonon-mode  eigen-frequencies  (from 
lowest  to  highest)  with  integers  (from  1  to  N  for  an  N-ion 
crystal)  and  using  a  linear  scaling  between  two  integers  to 
denote  detuning  located  between  the  corresponding  two 
modes.  For  instance,  in  this  notation  ft  =  2.75  means  the 
physical  detuning  q  =  002  +  0.75(co>3  —  u^)-  Each  phase 
is  characterized  by  a  spin  order  (denoted  with  a  binary 
string  where  0  and  1  correspond  to  t  and  4-  spin  respec¬ 
tively)  which  gives  one  of  the  ground  state  spin  configu¬ 
rations.  The  Ising  Hamiltonian  (1)  features  a  reflection 
symmetry  and  an  intrinsic  Z2  symmetry  with  respect  to 
a  global  spin  flip.  The  spin  order  breaks  the  Ising  symme¬ 
try,  so  each  phase  is  at  least  two- fold  degenerate.  If  the 
spin  order  also  breaks  the  reflection  symmetry,  the  corre¬ 
sponding  ground  state  is  4-fold  degenerate.  For  instance, 
for  the  phase  denoted  by  the  spin  order  01001,  the  four 
degenerate  ground  states  are  1 4 tilt)  5  |titti)>  |tiiti)> 
and  |ttttt)-  When  the  re-scaled  detuning  jl  crosses  an 
integer  (a  phonon  mode),  the  spin  order  changes  as  ex¬ 
pected,  but  this  is  not  a  conventional  phase  transition  as 
the  parameters  Jmn  change  discontinuously  in  the  Hamil¬ 
tonian  (1).  However,  when  jl  varies  within  two  adjacent 
integers,  all  the  parameters  Jmn  are  analytic  functions 
of  q,  yet  the  spin  order  can  still  change  abruptly,  sig¬ 
naling  a  phase  transition.  The  frequency  of  this  type  of 
inter- mode  phase  transition  increases  rapidly  with  the  ion 
number:  there  is  one  such  transition  for  a  three- ion  chain 
and  12  such  transitions  in  a  nine- ion  crystal.  Another  no¬ 
table  feature  from  Fig.  2(a)  is  that  there  is  typically  no 
phase  transition  when  jl  varies  from  an  even  mode  (2k) 
to  an  odd  mode  (2k  +  1,  k  =  0, 1,-  •  • ,  (N  —  1)  /2).  In 
such  regions,  the  spin  order  has  a  reflection  symmetry. 
This  suggests  that  a  spin  order  with  reflection  symmetry 
may  be  more  stable  in  energy  and  does  not  easily  yield 
to  other  spin  configurations.  This  observation  is  consis¬ 
tent  with  the  fact  that  for  an  even  number  of  ions,  there 
are  much  fewer  inter- mode  phase  transitions,  as  the  spin 
order  in  these  cases  has  a  reflection  symmetry. 

As  we  add  a  transverse  B  field  to  the  Hamiltonian, 
the  spins  will  gradually  become  polarized  along  the  x- 
direction  along  B.  In  Fig.  2(b),  we  plot  the  aver¬ 
age  polarization  )/N  as  a  function  of  the  field 

B  (in  the  unit  of  the  average  Jmn  defined  by  J  = 
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Figure  2:  (Color  online)  (a)  Ground- state  phase  diagram  at  B  =  0  characterized  by  the  corresponding  spin  orders  for  N  = 
3,5,  7,  9  ions.  The  transition  point  is  denoted  by  the  re-scaled  detuning  jl.  (b)  Average  polarization  for  N  =  7 

ions  at  finite  fields  B.  The  right  figure  is  a  closeup  near  the  critical  point. 


\Jmn\2/[N(N  -  1)])  for  N  =  7  ions  in  a  small 

region  of  the  detuning  jl.  We  find  that  the  system  is  eas¬ 
ily  polarized  if  it  lies  at  the  critical  point  between  two 
different  spin  orders  given  by  the  Ising  couplings.  But 
near  the  center  of  a  spin  phase,  the  spin  order  is  more  ro¬ 
bust  and  can  persist  under  a  finite  H,  eventually  yielding 
to  the  polarized  phase  as  B  increases  through  the  Ising- 
type  transition  (which  becomes  a  broad  crossover  for  this 
finite  system). 

With  B  =  0,  the  transition  between  different  spin  or¬ 
ders  is  sharp  as  it  is  characterized  by  a  level  crossing 
for  the  ground  state  of  the  Hamiltonian  (1).  When  we 
turn  on  a  finite  B  field,  the  system  shows  only  avoided 
level  crossings  in  its  ground  state,  and  we  expect  that 
the  sharp  phase  transition  at  B  =  0  to  be  replaced  by  a 
broad  crossover  for  this  small  system,  similar  to  the  Ising- 
type  of  transition  discussed  above.  Interestingly,  this  is 
not  always  the  case.  Even  at  a  finite  H,  some  transitions 
remain  very  sharp  in  this  small  system,  and  the  sharp¬ 
ness  increases  exponentially  with  the  number  of  particles. 
This  provides  an  unusual  finite-size  scaling  behavior,  dis¬ 
tinctively  different  than  the  polynomial  finite-size  scaling 
of  the  transition  width  that  one  typically  sees  in  conven¬ 
tional  phase  transitions  [7].  (For  instance,  the  conven¬ 
tional  Ising  transition  displays  a  linear  scaling  law  [7].) 
To  characterize  this  unusual  finite-size  scaling  behavior, 


we  look  at  a  particular  example:  for  N  ions,  there  is  a 
unique  spin  phase  transition  in  the  region  with  detuning 
N  —  2  <  jl  <  N  —  1  when  N  is  odd.  A  schematic  phase 
diagram  for  this  region  is  shown  in  Fig.  3(a).  At  B  =  0, 
we  have  a  ferromagnetic  phase  on  the  left  side  which  is 
doubly  degenerated  and  a  kink  phase  on  the  right  side 
which  is  4- fold  degenerate.  At  finite  T>,  these  two  spin 
orders  remain  robust  in  a  range  of  T>,  before  eventually 
yielding  to  a  polarized  phase  for  a  large  B  field  through 
a  crossover.  The  transition  between  the  ferromagnetic 
and  the  kink  phases  remain  sharp  for  a  finite  T>,  and 
this  sharpness  increases  rapidly  with  the  ion  number,  as 
one  can  see  from  Fig.  3(b)  and  3(c).  With  a  moderate 
increase  of  the  ion  number  from  5  to  9,  the  transition 
between  the  ferromagnetic  and  the  kink  phase  becomes 
almost  infinitely  sharp.  It  is  also  interesting  to  note  that 
the  transition  line  between  these  two  phases  has  a  slope 
with  the  B  axis,  so  one  can  cross  this  phase  transition  by 
tuning  either  the  detuning  jl  or  the  field  B. 

We  have  characterized  the  transition  width  between 
the  ferromagnetic  phase  and  the  kink  phase,  and  find  this 
transition  width  indeed  shrinks  exponentially  with  the 
number  of  ions.  To  see  this,  we  show  the  four  lowest  lev¬ 
els  across  this  phase  transition  when  we  tune  either  jl  or 
B  (see  Fig.  4a).  While  the  lowest  state  is  a  smooth  func¬ 
tion  of  jl  at  a  finite  T>,  the  second  and  third  states  feature 
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Figure  3:  (Color  online)  (a)  The  schematic  phase  diagrams  for  odd  numbers  of  ions  in  the  region  with  detuning  N— 2  <  jl  <  N—l. 
For  a  small  number  of  ions,  the  solid  line  represents  a  sharp  transition,  whereas  the  dashed  lines  represent  a  continuous  crossover 
to  the  polarized  state.  (b,c)  The  calculated  theoretical  phase  diagrams  for  (b)  N  =  5  and  (c)  N  =  9  ions.  Color  shows  the 
order  parameter  defined  by  Pfm  —  Pk ,  where  Pfm  and  Pr  are  the  projection  probabilities  of  the  ground  state  of  the  system  to 
the  Hilbert  subspace  with  the  ferromagnetic  and  the  kink  orders,  respectively  (the  basis- vectors  of  the  corresponding  subspaces 
are  indicated  in  Fig.  (a)). 


Figure  4:  (a)  The  structure  of  the  lowest  four  energy  levels 
around  the  sharp  ferromagnetic-kink  phase  transition.  The 
two  lowest  states  have  a  ferromagnetic  (kink)  order  on  the 
left  (right)  side.  The  transition  width  W  is  proportional  to 
the  energy  gap  A E,  and  AE  is  fitted  by  an  exponential,  (b) 
The  exponent  a  is  a  function  of  the  ion  number  N,  where  the 
solid  line  is  a  linear  fit  a  ~  ( N  —  1)  /2. 


a  level  crossing.  From  Fig.  4a,  the  transition  width  W 
is  proportional  to  the  energy  gap  AE.  We  calculate  the 
energy  gap  AE  as  a  function  of  the  ion  number  TV,  and 

find  it  is  well  fit  by  AE  oc  (^B/NJ^j  for  B  NJ ,  where 

the  exponent  a  has  a  roughly  linear  dependence  with  N 
as  shown  in  Fig.  4b  and  can  be  fit  by  a  ~  (TV  —  1)  /2. 
This  formula  shows  that  the  energy  gap  AE ,  and  thus 
also  the  transition  width  W,  shrinks  exponentially  with 


the  ion  number,  which  explains  why  the  transition  is  so 
sharp  for  a  small  system.  The  exponential  shrinking  of 
the  energy  gap  A E  can  be  intuitively  understood  as  fol¬ 
lows:  when  B  <C  TV  J  we  can  treat  the  term  B 
as  a  perturbation  in  the  Hamiltonian  (1).  For  each  ap¬ 
plication  of  B^2ncTni  we  can  only  flip  the  direction  of 
one  spin.  As  the  ferromagnetic  state  and  the  kink  state 
have  (. N  —  1)  /2  spins  taking  opposite  directions,  the  two 
states  need  to  be  connected  through  (N  —  l)  /2th  order 
perturbation,  and  thus  the  energy  gap  is  proportional  to 
(. B/Nj)(N~iy 2. 

In  summary,  we  have  shown  that  laser  induced  mag¬ 
netic  coupling  between  trapped  ions  realizes  a  frustrated 
Ising  spin  network  with  competing  long  range  interac¬ 
tions,  giving  rise  to  rich  phase  diagrams  for  the  ground 
state.  Some  of  the  phase  transitions  in  this  system  are 
characterized  by  a  unusual  finite  size  scaling,  where  the 
transition  width  scales  down  exponentially  with  the  num¬ 
ber  of  ions.  This  exponential  finite  size  scaling  lads  to 
sharp  phase  transitions  for  a  small  system  even  with  just 
a  few  ions,  as  one  can  realize  now  in  the  lab. 

This  work  was  supported  by  the  DARPA  OLE  pro¬ 
gram,  the  IARPA  MQCO  program,  the  AFOSR  and  the 
ARO  MURI  programs,  the  NSF  PIF  Program,  and  the 
NSF  Center  at  JQI. 


[1]  D.  Porras  and  J.  I.  Cirac,  Phys.  Rev.  Lett.  92,  207901 
(2004). 

[21  X.-L.  Deng,  D.  Porras,  and  J.  I.  Cirac,  Phys.  Rev.  A  72, 
063407  (2005). 

[3]  A.  Friedenauer  et  al.,  Nature  Physics  4,  757  (2008). 

[4]  K.  Kim  et  al.,  Phys.  Rev.  Lett.  103,  120502  (2009). 

[5]  K.  Kim  et  al.,  Nature  465,  590-593  (2010). 

[6]  E.  E.  Edwards  et  al.,  Phys.  Rev.  B  82,  060412  (2010). 

[7]  S.  Sachdev,  Quantum  Phase  Transitions,  Cambridge  Uni¬ 


versity  Press,  Cambridge  (1999);  P.  Ribeiro,  J.  Vidal,  and 
R.  Mosseri,  Phys.  Rev.  E  78,  021106  (2008). 

[8]  In  the  experiments  of  Refs.  [4-6],  the  lasers  drive  a  (in¬ 
dependent  force  (see  A.  Sorensen  and  K.  Molmer,  Phys. 
Rev.  A  62,  022311  (2000).)  Here,  we  use  the  standard 
basis  of  the  Ising  interaction  by  exchanging  cr*  and  cr* 
for  convenience,  but  this  simple  change  of  basis  does  not 
modify  the  essential  physics  in  this  system. 

[9]  G.-D.  Lin  et  al.,  Europhys.  Lett.  86,  60004  (2009). 


[10]  S.  -L.  Zhu,  C.  Monroe,  and  L.  -M.  Duan,  Phys.  Rev.  Lett.  Publishing  Company  (2005). 

97,  050505  (2006).  [12]  F.  Barahona,  J.  Phys.  A:  Math.  Gen.  15  3241  (982). 

[11]  H.  T.  Diep,  Frustrated  Spin  System,  World  Scientific 


